Field-dependent spin and heat conductivities of dimerized spin- 1/2 chains 
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We study the spin and heat conductivity of dimerized spin-1/2 chains in homogeneous magnetic 
fields at finite temperatures. At zero temperature, the model undergoes two field-induced quantum 
phase transitions from a dimerized, into a Luttinger, and finally into a fully polarized phase. We 
search for signatures of these transitions in the spin and heat conductivities. Using exact diagonal- 
ization, we calculate the Drude weights, the frequency dependence of the conductivities, and the 
corresponding integrated spectral weights. As a main result, we demonstrate that both the spin 
and heat conductivity are enhanced in the gapless phase and most notably at low frequencies. In 
the case of the thermal conductivity, however, the field-induced increase seen in the bare transport 
coefficients is suppressed by magnetothermal effects, caused by the coupling of the heat and spin 
current in finite magnetic fields. Our results complement recent magnetic transport experiments on 
spin ladder materials with sufficiently small exchange couplings allowing access to the field-induced 
transitions. 
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I. INTRODUCTION 

While low-dimensional quantum magnets have been in- 
tensely studied for the past decade, the transport proper- 
ties still pose viable challenges for experimental and theo- 
retical physicists ii"— On the theoretical side, the ground- 
state properties of low-dimensional spin systems are very 
well studied by means of powerful techniques such as the 
Bethe ansatz, bosonization, the density matrix renormal- 
ization group (DMRG) method, quantum Monte Carlo 
(QMC), and exact diagonalization (see Refs. la-Q for re- 
views). Finite-temperature transport properties, though, 
remain an exciting and active field of research. One 
of the best established results is the ballistic thermal 
transport in the integrable XXZ spin-1/2 chainj^"— In 
general, however, transport properties are not that eas- 
ily obtained, especially far from equilibrium or in non- 
integrable models, but even spin transport in exactly 
solvable models is still an active field of research with 
open pending questions j^i"— In the non-equilibrium case, 
one has to resort to numerical simulations of either sys- 
tems coupled to external bathai^Ti^i or closed systems 
prepared in a state out of equilibriumj ^^i^'^ 

For non-integrable systems, finite-temperature trans- 
port is expected to be diffusivC )^^'^^ which is consistent 
with the numerical observation of a vanishing Drude 
weight for both heat and spin transport j ^^i^^" — While 
this result has been accepted for massive phases in 
the high-temperature regime, the situation for massless 
phases of non-integrable systems at low temperatures and 
in the vicinity of integrable modelsiii^Sr— is less clear, in 
the sense that numerical studies yield large ballistic con- 
tributions to transport coefficients in massless phases and 
on finite systems i^ '^^i"^^ 

Theoretical research into the transport properties 
of low-dimensional magnetsi*^ has been strongly moti- 
vated by exciting experimental results on materials with 
low-dimensional electronic structures^i^ The large ther- 



mal conductivities found in the spin ladder materials 
(Sr,Ca,La)i4Cu2404i establish a link between the ther- 
mal conductivity and magnetic excitationsi^"— There 
is also a variety of spin chain materials with magnetic 
contributions to the heat conductivity, most notably 
SraCuOg, SrCuOs (Ref.!!!-©, and CaCuaOg (Ref.Ell). 
Recent successes in sample preparation have resulted in 
very clean samples of SrCu02, exhibiting the largest 
thermal conductivity so far observed in low-dimensional 
quantum magnets4^ This has been interpreted as ex- 
perimental evidence for ballistic heat transport in clean 
Heisenberg chains, where phonons are the main source of 
external scattering4£ 

The magnetic field dependence of the thermal conduc- 
tivity has been the case of interest in several experi- 
mental studies^i^i^ yet in most of the known mate- 
rials, exchange couplings are orders of magnitude larger 
than the magnetic fields available in a laboratory. Only 
recently, the quasi-one-dimensional organic compound 
(C5Hi2N)2CuBr4 (Refs. lillisl) has gained attention in 
this context. It has exchange couplings small enough to 
allow experimental access to field-driven quantum phase 
transitions at low temperatures up to the saturation 
field. In these experiments, the phase diagram with re- 
spect to temperature and magnetic field has been ex- 
plored by a large variety of experimental probes, estab- 
lishing the presence of a field-induced gapless phase at 
low temperatures4^"— Recent measurements of the ther- 
mal conductivity of these compounds in external mag- 
netic fields have been interpreted in terms of the absence 
of spin mediated heat transport)^ 

The field-dependent thermal transport in the XXZ 
chain has previously been addressed with several theo- 
retical approaches)^""— emphasizing the role of magne- 
tothermal corrections to the thermal conductivity due 
to the coupling of the spin and the heat current, simi- 
lar to the Seebeck effect i^ The possibility of controlling 
the heat transport in spin chains by varying a magnetic 



field has been addressed in Ref. \5Q. The zero-field trans- 
port properties of the dimerized chain have been studied 
in Refs. IMllI while the field-dependence of the thermal 
Drude weight of noninteracting, dimerized XX chains has 
been discussed in Ref. [53- 

Our present goal is to understand the dependence of 
the spin and heat conductivities of dimerized spin chains 
on external magnetic fields. Within linear response the- 
ory, the frequency-dependent conductivity has two con- 
tributions, a delta peak at zero frequency whose weight 
is the Drude weight and a regular part: 



(t[k]{uj) = DsHi,]S{l^) + Cr[K]rog(w) • 



(1) 



Using exact diagonalization to evaluate Kubo 
formula o^^'^^ we compute these quantities for a value 
of the spin gap that is comparable to the one found in 
(C5Hi2N)2CuBr4 (Refs.lHl^). 

As a main result, we find an increased weight in the 
low-frequency regime of both <7{lo) and k{uj) in the field- 
induced phase. This allows for a direct interpretation 
in terms of transport channels opened in the vicinity 
of the Fermi points of the corresponding Luttinger liq- 
uid. These channels lead to a strong enhancement of the 
transport coefficients for spin and heat conduction. On 
the finite systems that we have access to with exact di- 
agonalization, the main contribution to the increase of 
the conductivities at low frequencies and low but finite 
temperatures stems from the Drude weight. In the case 
of the spin conductivity this increase is roughly an order 
of magnitude larger than the one observed in the regular 
part, whereas for the thermal conductivity, the picture 
is more involved. At zero magnetic field and low tem- 
peratures, the bare thermal Drude weight (without any 
magnetothermal corrections) dominates the regular part, 
yet when increasing the field, the regular part increases 
more strongly than the thermal Drude weight. Finally, 
taking into account the magnetothermal corrections, the 
increase of the thermal weight with increasing field be- 
comes a decrease, as expected from the results for the 
XXZ chain,^ 

We further study the dependence of the transport coef- 
ficients on the strength of the dimerization, varying it be- 
tween the limits of uncoupled dimers and the Heisenberg 
chain. We expect that our results are generic for dimer- 
ized quasi one-dimensional systems, while the obvious 
advantage of working with the dimerized chain is that 
with exact diagonalization, we can reach longer chains 
than in the case of a ladder. We shall stress that our 
work is concerned with the intrinsic transport properties 
of dimerized systems, whereas for a complete description 
of the experimental results, phonons may play an impor- 
tant role, as has been emphasized in Refs. [59l - l6ll . 

The paper is organized as follows: First, we introduce 
the model and briefly review the ground-state proper- 
ties. Second, we proceed by summarizing the necessary 
framework to compute transport coefficients and conduc- 
tivities within linear response theory. Section IV presents 
the results. We study the Drude weights in Sec. IIV Al 
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FIG. 1: (color online) (ft, A)-plane of the phase diagram of the 
dimerized spin chain at temperature T = (see Refs. KJamS ). 
The upper (red dashed) line is the saturation field. The lower 
(solid black) line corresponds to the spin gap.— "^ The critical 
fields ftci and hc2 are obtained from exact diagonalization of 
chains of length L = 18 (squares). The blue diamonds are 
DMRG data, extrapolated in the inverse system size. 



the frequency-dependent conductivities in Sec. IIV Bl and 
the spectral weights in Sec. IIV CI all as a function of 
the magnetic field and temperature. The influence of 
the strength of the dimerization is discussed in Sec. IIVDI 
and Sec. IIV El covers magnetothermal effects. Finally, we 
summarize our findings in Sec. |V| The dependence of 
the Drude weights at finite fields on the system size is 
presented in the Appendix. 



II. MODEL 

We study the spin and thermal conductivity of dimer- 
ized spin chains in a homogeneous magnetic field. The 
model Hamiltonian for a chain of length L is given by 



H 



^/ij = j'^[XiSi ■ Si- 



hSf 



(2) 



where the dimerization is introduced via: Ai = A for 
i even and A, = 1 for i odd, S = {S'',Sy , S""), Sf and 
/i = x,y,z are the components of a spin-1/2 operator 
acting on site i and h denotes the magnetic field. We use 
periodic boundary conditions. The exchange coupling 
J sets the global energy scale and the model includes 
the Heisenberg spin chain (A = 1) and uncoupled dimers 
(A = 0) as limiting cases. 

The magnetic phase diagram^^i^ shown in Fig. [U is 
very similar to the one of the two-leg ladder systemi^ 
The upper (red dashed) line is the saturation field. The 
lower (solid blue) line corresponds to the spin gap which 
has been intensely studied by analytical and numerical 
meansi^"— Since the Heisenberg chain is gapless a care- 



ful finite-size scaling for tlie lower critical field is nec- 
essary. The blue diamonds are DMRG data for open 
boundary conditions, extrapolated in inverse system size, 
which we compare to exact diagonalization results ob- 
tained with L = 18 sites and periodic boundary condi- 
tions (squares). This illustrates that for A < 0.8, the 
phase boundaries exhibit very small finite-size effects; 
thus, for the system sizes that we shall use in our ex- 
act diagonalization analysis, we are already very close to 
the bulk value for the gap. The Luttinger liquid phase 
here is similar to the Luttinger liquid phase in the XXZ 
model where the transition from the gapped to the gap- 
less phase is of the commensurate-incommensurate type 
(see, e.g., Refs.lHandllS). 



III. TRANSPORT COEFFICIENTS AND 
CONDUCTIVITIES 

Here we summarize the central equations for magnetic 
transport in the linear response regime. The expectation 
values of the spin and thermal currents, ji and j2, are 
given by^ 



{Ji 



) ' / ^ ^Iriijn 



(3) 



where /i = Vh and /2 = — VT refer to the magnetic 
field and temperature gradients. Lim is the conductivity 
matrix, j'l and J2 can be expressed via the spin and 
energy currents js and jth by 



where 



jl=js, J2^jth~hjs 
L-l 

js[th] = i^[hi-i,di] . 
1=1 



(4) 



(5) 



hi denotes the local energy densities defined by Eq. ([2]) 
at zero magnetic field, i.e., di = hi = XiS^S^_^i for heat 
transport, while di = S[ for spin transport. The conduc- 
tivities satisfy 



0';,n(w) = ReL/„i (w) = DirnS{uj) + CTrcg,;7n('^) (6) 



with 



A, 



7r/3' 



T + 1 



ZL 



J2 e-^''"{n\ji\o){o\j^\n), (7) 



E„=Ea 

where r = 0(1) for m = 1(2) and 
_ TT^'' 1 - e-'3" 



E 



,-/9-E„ 



ZL uj 

X (n|ji|o)(o|j„»(5(cj-AS). (8) 

/3 = 1/T is the inverse temperature, \n) and E^ are the 
eigenstates and energies oi H, Z = ^^ q-PE„ (denotes 



the partition function, and AE 
define 



Eo — En- Next we 



K (w) 

crth,s (w) 



CTll (W) 

0-21 {i^)\ 



(9) 



j2->Jth 



and the Drudc weights Dg. Dth, and Dth.s of cr(a;), ^(a;), 
and (Tth,s('^)j respectively. The latter set of quantities 
corresponds to a choice for the currents alternative to 
Eq. ((4]), namely ji = js and J2 = itti- Note that 



TD 



12(s,th) 



D2I 



(th,s)- 



As in Refs. Ulllliii we will 
evaluate a (w), k (uj), and crth,s (w) using exact diagonal- 
ization. From these quantities all transport coefficients 
o'lmi^) at finite fields can be obtained using Eqs. (UJ, ([71), 
and ([H]) . The Drude weight Kth for purely thermal trans- 
port accounts for the situation of no spin current, i.e., at 
(js) = 0. It is obtained analogously to the Seebeck effect 
as 



K< 



th 



D 



th 



D2 



(10) 



or equivalently, using Eqs. ^ and ([7]): Kth — D22 — 
D2i/{TDii). We will study the field dependence of the 
magnetothermal coupling in Sec. IIVEI 

For the spin Drude weight, an alternative expression 
fully equivalent to Eq. ([7]) exists^! 



^^-ZL^^-^^ 






-;3£;„ K47s|n)f 



AE 



(11) 



where (T) denotes the kinetic energy. This expression 
will be used to evaluate Dg in this paper since it gives 
the correct contribution to the optical sum rule Eq. (|13p 
on finite systems at low temperatures (see the discussion 
below) ,^ 

Finally, we define the integrated spectral weights 

a[K]{uj)duJ = Ds[thl + 2 / Cr[K]rcg(t^) , 

(12) 
and I^ui-A = /s[th](oo). For the spin conductivity one ob- 
tains the optical sum rule^^ 

I?:-- 



(j{uj)duj = -(-T). 

1j 



(13) 



The r.h.s. of the corresponding sum rule for thermal 
transporli^^ at finite temperatures depends on the model 
and the choice for the local energy density. 



IV. RESULTS 

We present our results for the transport coefficients 
of dimerized spin chains, starting with the dimeriza- 
tion strength A = 0.5, yielding a gap close to the 
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FIG. 2: (color online) Drude weight for spin (a) and thermal (b) transport at A = 0.5 and different temperatures for chains 
of length L = 18. At low temperatures (T = 0.25J, solid black lines) we find a significant increase of both quantities in the 
field-induced phase. At larger temperatures, the peak becomes less pronounced and is located at higher magnetic fields. The 
black dotted vertical lines are guides to the eye which mark the critical fields hc^ = 0.66J and hc2 ~ 1.4J. 




FIG. 3: (Color online) Regular part of the dynamical conductivities cr(aj), k(cj) as a function of frequency for different magnetic 
fields and temperatures, (a) shows the spin conductivity at low temperatures (T = 0.25J) and zero magnetic field (dashed 
black line) as well as h = J (solid red line) for a chain with L = 18 spins and a dimerization of A = 0.5. In the presence of 
a field h = J, crreg(w) increases in the gap, but strongly decreases at higher frequencies. The inset (b) shows the same at a 
higher temperature (T = J) where the effect of h is less strong, (c) and (d) show Krcg(^) of the same system, (c) For low 
temperatures, the magnetic field gives rise to an increase of Krcg(^) over the whole spectrum (solid red line). The inset (d) 
shows Krog(i^) at higher temperatures (T = J), where Krcg(i^) is a smooth curve and barely influenced by the fleld. 



one of the experimental system. We separately ana- 
lyze the Drude weights (Sec. lIVA| and the conductivities 
(Sec. IIVB[) and then combine them to get the spectral 
weight (Sec. lIVCp . all as a function of temperature and 
magnetic field. As a result, we find that the dominant ef- 
fect of applying a magnetic field for this fixed value of A is 
an enhancement of the low-frequency spin as well as heat 
conductivity. In Sec. IIVDI we analyze this enhancement 
at different strengths of the dimerization. Section IIVEI 
discusses the magnetothermal coupling. 



A. Drude weights as a function of T and h 



We start focusing on A — 0.5, resulting in a value 
of the spin gap A = 0.66 in units of J. One of the 
two contributions to the integrated spectral weights is 
the Drude weight in Eq. ([T]). Figure [5] shows the spin 
Drude weight Dg and the thermal Drude weight Dth as 
a function of the magnetic field for T/J = 0.25,0.5,1 
and chains consisting of L = 18 spins. At low tempera- 
tures (T = 0.25 J, solid black lines), we find a strong en- 
hancement of both quantities between the critical fields 
/ici = 0.66J and hc2 — 1.48J. Increasing the tempera- 
ture weakens this feature significantly, especially for ther- 
mal transport. This feature is addressed quantitatively 




FIG. 4: (Color online) (a) Spectral weight /s(aj) for spin transport integrated up to cj = 0.5 J as a function of magnetic field h 
for L = 18 spins and at different temperatures, (b) Spectral weight Ith{^) for heat transport integrated up to w = 0.5 J. At 
low temperatures, there is a large peak in both quantities indicating the enhancement of transport in the field-induced gapless 
phase. As temperature increases this effect vanishes up to the point where the magnetic field causes the decrease of spectral 
weight in the gap. For heat transport, the magnetic field has almost no effect at large temperatures (T = J, 2J). (c), (d) When 
extending the upper bound on the frequency oj to infinity, the effect of the magnetic field on spin transport is hidden under 
the large contributions of the regular part above the gap (c), while I^^^ behaves similar to Ith{^ = J/'2) (d). The inset in each 
panel illustrates the small dependence on the system size at low temperatures. 



in the discussion of the spectral weights (Sec. IIVC|) . It 
is also worth noting that the maximum in the Drude 
weights i'sfth] = DsUY^]{h) moves to higher fields as the 
temperature increases. At low temperatures, the ques- 
tion arises whether there are any features in Dg = Ds{h) 
at the critical field, i.e., whether the finite systems we in- 
vestigate exhibit remnants of the quantum critical behav- 
ior to be expected at hc^ . We shall comment on this point 
in Sec. IIVDI where we will discuss the A-dependence of 
the Drude weights. 



B. Dynamical conductivities 

The second contribution to the spectral weight is the 
regular part [Eq. (|5])] of the dynamical conductivity. By 
inspecting the spectral representations, one realizes that 
the field can only enter via the Boltzmann factors, thus a 
magnetic field can only enhance or reduce existing weight 
but does not change the position of the poles cu = AE. 



Note that at zero magnetic field, the weight in the gap 
is suppressed at low temperatures as the system ap- 
proaches its dimerized ground-state and can be tuned 
to finite values by increasing either the magnetic field or 
temperature 1^ 

Our numerical results are depicted in Fig. [31 At zero 
magnetic field and low temperatures, there is almost 
no weight in the regular part of cr(w) [dashed line in 
Fig. [3l^a)] below a value of to approximately correspond- 
ing to /ici (compare Ref. for the case of A = 0.1 and 
h = 0) beyond which the dominant peak is located. Turn- 
ing on a magnetic field larger than hc^ influences the 
curve drastically. The major portion of the weight still 
lies above the gap but is much smaller and without sig- 
nificant peaks [solid line in Fig.[21[a)]. In addition to the 
overall reduction we find spectral weight at very small 
frequencies which is not the case at zero field. Going 
to higher temperatures [Fig. ^h)] results in a smooth 
curve due to thermal excitations, while the influence of 
the magnetic field is much weaker and does only change 




FIG. 5: (color online) Relative enhancement A7a[th] of the 
spectral weights as a function of temperature (L = 14). The 
curves for thermal transport (solid black line) and spin trans- 
port (dashed red line) almost coincide and decay slightly 
slower than exponentially (A/s[th] ~ exp(— cT), green dash- 
dotted line) as the temperature increases. 



the numerical values without modifying the structure. 

In the case of thermal transport [Fig. El^c)], the basic 
structure is different and the role of the magnetic field 
at low temperatures is even more important. Without a 
field (dashed line) a significant amount of spectral weight 
is found around a; = while another important contribu- 
tion is at the lower edge of the one-triplet band, similar 
to the spin conductivity. Switching on a magnetic field 
larger than h^ (solid line) enhances the heat conductiv- 
ity strongly, especially at low frequencies. Increasing the 
temperature up to T = J smoothens the curve while com- 
pletely suppressing the infiuence of the magnetic field. 

Concentrating on the regular parts alone, one would 
conclude that both spin and heat transport are not infiu- 
enced by the magnetic field at high temperatures. At low 
temperatures the weight in o-i.cg(a;) is strongly reduced by 
the magnetic field yet comes along with an increase in the 
weight at low frequencies, which is the effect we are inter- 
ested in. By contrast, Ki.eg(w) is enhanced by the field at 
low temperatures at all frequencies, with most of this in- 
crease originating from the low-frequency range to < h^^ . 
At higher temperatures, both effects are suppressed. To 
understand the interplay of the Drude weights and the 
regular parts one has to study the spectral weight. 



C. Field-dependence of the spectral weight 

Now that we have shown how the transport enhance- 
ment takes place via increasing weight at low frequencies 
we study the integrated spectral weight /s[th](a;). Fig- 
ure |4] shows our results for Is{uj = 0.5 J) and Ith{^ = 
0.5 J) as well as the total spectral weights I^^m- We have 
chosen the upper bound of integration to be a; = 0.5 J 
to measure how the weight in the gap of the system be- 



haves as a function of temperature and field. Our find- 
ings in the case of low temperature and A = 0.5 are the 
following: For both spin and heat transport, we find a 
strong enhancement of transport in the field- induced gap- 
less phase at low temperatures T < hc^ . As temperature 
increases, the field dependence becomes less pronounced. 
From Fig. HI it is obvious that the total weight in the 
spin conductivity, i.e., I^ oc (— T) typically decreases as 
a function of increasing field, due to the increasing weight 
of strongly polarized states. Therefore, the enhanced low- 
frequency weight seen in both the Drude weight and the 
regular part has to be accompanied by a decrease of spec- 
tral weight at higher frequencies in order to satisfy the 
optical sum rule Eq. (|13p , consistent with our discussion 
ofFig.[3i;a). 

In the case of thermal transport, the infiuence of tem- 
perature on I^^ is much more drastic: the maximum in 
/°jj seen at T = 0.25 J is quickly washed out at tempera- 
tures T > hc^. By contrast, for spin transport and in /^, 
the interesting features are hidden at all temperatures 
under the large field-induced decrease of the regular part 
as w approaches the gap [Fig. [^c)], whereas in the case 
of thermal transport, the low- frequency behavior domi- 
nates the field dependence of I^^ at T < hc^ [Fig. IHJd)] . 
For T > 0.25J, the finite-size effects of the integrated 
weights become very small, which is shown in the insets 
of Fig. |4l illustrating that this quantity is very robust. 
In other words, the maximum of Is[^^{Ld ~ 0.5 J) in the 
window /ici < h < hc^ is stable against a variation in L. 

The field-dependent enhancement of /sfthll"^ = 0.5 J) 
in Fig. H (a) and (b) (T = 0.25J) is weU described by 



s[th] 



h^ for h < hci . Furthermore we present 



_ max;t(/s[th](^ = 0-5 J, h)) - /s[th](^ = 0-5 J, fe = 0) 
''*''' ~ /s[th](w = 0.5J,/i = 0) 

in Fig. [5] as a function of temperature for L ~ lA spins. 
The data are already converged with respect to the sys- 
tem size for the temperatures shown. While the field- 
dependence of Ith{i^ = 0.5J) and /s(w = 0.5J) is only 
qualitatively similar, A/th (black solid line) and A/g (red 
dashed line) almost coincide. The green dash-dotted line 
is an exponential decay fitted to AJth to illustrate that 
the decay is slightly slower than exponential for spin as 
well as for heat transport. 

Regarding the temperature dependence of the thermal 
spectral weight at lower temperatures (T < 0.25J) we 
expect an exponential suppression at low temperatures 



in gapped phases 



11.57.74 



{h < he 



h > hr 



For 



the Luttinger Liquid phase {h^ < h < hc^), we expect 
-^s[th] — -^s[th] (T) ^ T based on previous work on the tem- 
perature dependence of the thermal drude weights of the 
XXZ spin- 1/2 chain in a magnetic field.— At zero mag- 
netic field we indeed find an exponentially suppressed 
weight converged with respect to the systems size for 
T < 0.25J, but in the gapless phase {h = J) the finite- 
size effects do not allow for a definite conclusion on the 
temperature dependence. 
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FIG. 6: (Color online) Drude weights normalized by the total spectral weight /^thi as a function of magnetic field, (a) In the 
case of spin transport, the Drude weight is a small contribution at zero field, but increases with the magnetic field up to a 
saturation value that strongly depends on temperature, (b) For heat transport the Drude weight is the dominant contribution 
even at zero field and low temperatures. At low temperatures the normalized Drude weight has a minimum as the field increases. 
At h ~ 2 J, it reaches its final value. Such a minimum is not visible at T = J. The black dotted vertical lines are guides to the 
eye which mark the critical fields hci = 0.66J and /ic2 = 1.4J. 
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FIG. 7: (Color online) Spin Drude weight as a function of magnetic field h for a system of length L = 16 and different values 
of A. (a) At low temperatures (T = 0.25J) the spin Drude weights show a significant peak between the critical fields (compare 
Fig. [T]at T = 0). (b) For larger temperatures (T = J) the peak is less pronounced and moves to larger fields. In both cases 
increasing A gives rise to a strong increase of the numerical values but the position of the peak moves only slightly towards 
smaller fields. The qualitative field dependence is almost independent of A. For all values of A the peak is located at marginally 
smaller fields than for the Heisenberg chain [long dashed lines in (a)]. 



To clarify the role of the Drude weight on finite sys- 
tems in this context, it is plotted in Fig. IH] as a fraction 
of the corresponding total spectral weight. In the case of 
spin transport [Fig. [5] (a)] , the Drude weight contributes 
very little to the total weight but its significance increases 
monotonously with the field. Ds/I^ saturates at a finite 
value above the saturation field. These observations hold 
for all temperatures studied {T/J = 0.25, 0.5, 1), only the 
value at saturation and the field necessary for reaching 
it depends on temperature. For instance, at T — 0.25 J 
(solid black line) the Drude weight accounts for 60% of 
the total weight above h = 2J. At T = 0.5J (dashed 
red line), Ds/I^ saturates at h = 3.5 J at a value of 48%. 
However, Fig. [5] should not be interpreted in terms of a 



large absolute value of the Drude weight above satura- 
tion in the thermodynamic limit, neither at zero or finite 
temperatures. At zero temperature, for N -^ oo and in 
the two gapped phases h < hc^ and h > hc2, we expect 
the spin Drude weight to vanish, according to Kohn's 
reasoningjii As far as the gapless phase is concerned, the 
zero temperature Drude weight should consequently be 
finite in the thermodynamic limit, yet the contribution 
of the Drude weight relative to the regular conductivities 
remains an open issue for low temperatures. 

In the case of heat transport [Fig. [6jb)], the picture 
has more facets. For low temperatures (T ~ 0.25 J, solid 
black line) Dth/Iti^ decreases in the field-induced phase 
as there is a huge increase in the regular part of the con- 



ductivity [compare Fig. |3l^c)] . This results in a minimum 
in D^i^/Itii at /i w J, roughly where the field-dependent 
thermal Drude weight has its maximum. At ft, w 2 J, the 
Drude weight is restored as the dominant contribution 
at Uth/^^th ~ 0.9 at low T. This effect is weakened by 
temperature [compare Fig.[3l^d)], eventually resulting in a 
monotonous increase of Dth/Ith with hatT—J. Our re- 
sults clearly suggest that the Drude weight is a significant 
contribution in the field-induced gapless phase on finite 
systems. Several comments are in order. First, a finite 
Drude weight at T > in a non-integrable system would 
be surprising, while on the other hand, this observation 
is consistent with the fact that thermodynamic proper- 
ties in the phase are well described by an effective XX Z 
uiodelr^^ which is believed to have ballistic transport 
properties in its gapless zero-field phasci^ii^ Moreover, in 
the ferromagnetic phase h > hc^ , the physics is expected 
to be well approximated by a weakly-interacting gas of 
magnons which in turn renders the Drude weight large 
in relation to the total weight. We stress that similar 
to Ds, this is a statement about finite systems and tem- 
peratures, and much larger system sizes (see Appendix) 
would be required to extrapolate the ratio Z^th/^^h ^'^ ^^'^ 
thermodynamic limit. We can therefore not pursue this 
question here. Note though, that in the high-temperature 
limit (/3 = 0) the magnetic field dependence drops out. 
In this limit, exact diagonalization yields a systematic 
decrease of the Drude weight with the system sizeiii 

To summarize the discussion of the behavior at A = 
0.5, in the low-temperature regime, which is our main 
case of interest, the field-driven enhancement is visible 
in the conductivities as well as in the Drude weights 
and therefore also in the spectral weight, but the Drude 
weights contribute the most. Further, we can distinguish 
a field- and a temperature-dominated regime: qualita- 
tively, at T < ftcj , a variation of the magnetic field infiu- 
ences the conductivities, whereas for T > hc^ , the mag- 
netic field has little effect on both the structure and the 
weight in the conductivities. Figure |4] contains the main 
result of our work: a field-induced increase in both trans- 
port coefficients at low frequencies uj < /ici ■ 



D. Dependence on A 

Next we turn to the discussion of the influence of the 
strength of dimerization A. Since we have presented evi- 
dence that the Drude weights are the dominant contribu- 
tion in the field-induced phase, we concentrate on these 
quantities, expecting them to reflect the main qualita- 
tive behavior. Figure [7] shows the spin Drude weight 
as a function of the external magnetic field at different 
temperatures and strengths of dimerization. In all cases, 
the Drude weight for spin transport exhibits a maximum 
at intermediate fields. We expect that in the thermo- 
dynamic limit and at sufficiently low temperatures, the 
transition into the field-induced Luttinger phase should 
lead to signatures in the Drude weights at hc^ and hc^- 



Yet, for the system sizes we can study, and thus the ac- 
cessible T, the location of the inflection points of Ds do 
not exhibit a clear correlation with critical fields ftci and 
hc2 of Fig. [TJ This remains to be analyzed in the future. 
As temperature increases [compare Fig. [7]^a) and 
Fig. El^b)], the maximum in the Drude weights moves 
from between the critical fields to higher values of h while 
the absolute value of Dg decreases. Sending A — > de- 
creases the values of D^ for all fields, while the field- 
driven increase is not affected qualitatively. This is true 
for small as well as high temperatures. Note that increas- 
ing the temperature suppresses the spin Drude weight 
much more severely than altering A. In the case of ther- 
mal transport [Fig. [5]^a)] we find a stronger quantita- 
tive dependence on A while the fact that altering A does 
not change the field dependence remains qualitatively 
correct. The former is not surprising, since uncoupled 
dimers cannot carry any heat current while the Heisen- 
berg chain is known for its ballistic heat transport.—"— 



E. Magnetothermal couplings 

The last effect to investigate is the infiuence of the 
magnetothermal coupling Eq. pU)) . In the case of the 
XX Z chain, it has been shown that this correction sup- 
presses the thermal conduct ivityi^ The results for the 
field-dependence of the thermal Drude weight iCth at low 
temperatures {T ~ 0.25 J) arc shown in Fig. [HUb). The 
corrected Drude weight Kth [Eq. p^ ] decreases as the 
magnetic field increases, which is the systematic behav- 
ior in the whole field-induced phase. Shifting A towards 
A = lowers the overall values while only slightly affect- 
ing the field dependence. Figure [His a comparison of the 
three Drude weights associated with thermal transport, 
Dth,D22 and A'th, where 



Do 



D 



th 



2(3hD, 



th,i 



fih^Ds . 



(14) 



As expected from Eqs. pOj) and (|T4|) . all three quanti- 
ties coincide at zero magnetic field. While D22 > £'th 
both quantities exhibit a similar field dependence , i.e., 
it features a maximum in the field-induced gapless phase 
and a smooth decay to zero as h increases beyond satu- 
ration. The magnetothermal correction changes this be- 
havior completely, yielding a monotonous decrease as the 
field increases. 



V. CONCLUSIONS AND SUMMARY 

In this work we studied spin and heat transport in 
dimerized spin— 1/2 chains in a magnetic field, in de- 
pendence of temperature T, magnetic field h and the 
strength of dimerization A. Focusing on the field depen- 
dence in the case of A = 0.5 we found that the transport 
coefficients at low frequencies for both, spin and heat 
transport, are strongly enhanced in the field- induced gap- 



_ (a)T=0.25J 




y 


N 




1 






X=0.2 


- 


- 


/ 




\ 




-- ?i=0.5 


- 


/ 
/ 






\ 




■-■ ?L=1 


- 






_ 


/ 






'\ 


- 


/ 






\ 




/ y"' 






^ \ 




/ / 






\ 




~ y- 


-■ 


-- 


\ 

1 2L^ — 





0.8 



0.6 



^ 



0.4- 



0.2 



N 






_ •. (b) T=0.25J 

\ 

\ 


X=0.2 

-- X=0.5 
■-■ X=0.8 
■-■ X=l 


- 


\ 


, . J 1 


- 



2 

h/J 



2 
h/J 



FIG. 8: (Color online) (a) Thermal Drude weight Dth as a function of magnetic field h for a system of length L = 16 and 
different values of A at low temperatures (T = 0.25 J). While the thermal Drude weight decreases fast as A decreases, the field 
dependence is not affected, (b) Corrected thermal Drude weight A'th- We find that the magnetothermal coupling prohibits the 
strong increase of the thermal Drude weight in the field-induced phase. 



less phase. For spin transport at low but finite tempera- 
tures the Drude weight becomes the dominant contribu- 
tion as the field increases. We stress that this an obser- 
vation for finite systems. While one may expect that the 
Drude weight remains relevant in the field-induced phase 
in the zero temperature,— thermodynamic limit, the the- 
oretically interesting question of a finite Drude weight at 
low but finite temperatures is beyond the scope of this 
work. 

In the case of heat transport, the emerging picture is 
more involved. Increasing the magnetic field up to satu- 
ration, we find that the regular part of the conductivity 
is vastly enhanced at all frequencies. Although the ther- 
mal Drude weight has a field-dependence similar to its 
spin counterpart, this leads to a decrease of the relative 
contribution in the spectral weight at low temperatures. 
However, the thermal Drude weight remains the domi- 
nant contribution on the finite systems studied here. In 
both transport channels, a temperature larger than the 
spin gap severely weakens all features related to the mag- 
netic field. 

Our main finding is the field-induced increase of both 
transport coefficients in the low-frequency window, which 
is robust even on the finite chains accessible to exact di- 
agonalization. While most of our work was concerned 
with parameters that mimic the energy scales typical for 
(C5Hi2N)2CuBr4 [Ref. IH-liil, we further studied the de- 
pendence on the strength of dimerization. There we ob- 
serve that altering A changes the overall values while the 
field dependence remains qualitatively the same. Thus 
we conclude that our observation of field-enhanced spin 
and heat transport is valid in dimerized spin chains. We 
also expect it to be more generally valid in other dimer- 
ized quasi one-dimensional spin systems, such as two-leg 
spin ladders. 

Finally, the calculation of the magnetothermal cou- 
pling at low temperatures unveils a very interesting cf- 



3- 



t^ 



. 2- 






1 1 1 
T=0.25J, >i=0.5 






— D,, 


- 


- 




— D,H 




: A 


- K 


- 




y>----\ 




- 


~ -^ 1 ^ ^i-----^ 


-~-.4-_ 





2 
h/J 



FIG. 9: (color online) Comparison of the various thermal 
Drude weights Ah, D22 [Eq. and [Eq. ^], Kth [Eq. ([TOll] 
for thermal transport at A = 0.5 and T — 0.25 J for a chain of 
length L = 16. While the three curves coincide ai h — 0, the 
field dependence exhibits huge differences. On the one hand, 
D22 given by Eq. (|14p (solid black line) shows a much stronger 
enhancement than Dth (dashed red line). On the other hand, 
including the magnetothermal correction (dot-dashed green 
line) leads to a monotonous decrease of Kth as the magnetic 
field increases. 



feet, similar to the XXZ chaini^ The corrected ther- 
mal Drude weight shows no increase in the field induced 
phase, instead, it decreases. 

In conclusion, we complemented several experimen- 
tal and theoretical studies which characterized the field- 
induced gapless phase )^"— by emphasizing here that 
clear fingerprints of this transition are present in current- 
current correlation functions and should thus, in princi- 
ple, manifest themselves in transport measurements. The 
spin conductivity is inaccessible at the moment, but 
can be extracted from quantities measured in NMR 
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FIG. 10: (color online) Finite-size effects in the field-dependent spin Drude weight (a) and thermal Drude weight (b) at two 
different temperatures {T = 0.25 J, J) for the system sizes L = 14, 16, 18. 
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Appendix A: Finite size scaling of the Drude w^eights 

The finite-size scaling of the Drude weights in non- 
integrablc spin chains has been previously studied at zero 



field and in the limit of infinite temperatureii without 
evidence of a finite Drude weight in the thermodynamic 
limit for non-integrable systems. Note that the field de- 
pendence drops out at /3 = 0. These studies include the 
dimcrizcd chain. We present our data for finite magnetic 
fields in Fig. [TU] Due to the interplay between the regular 
contribution and the Drude weight in the conductivities, 
the dependence on the system size is non-monotonic at 
finite magnetic fields and for the system sizes accessible, 
especially at low temperatures (T = 0.25 J). Therefore, 
the accessible system sizes are too small to gain any qual- 
itative insight into -Ds[th] beyond Ref. [Ill- The integrated 
spectral weight (Fig. |3]) is a robust quantity in this con- 
text. 
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